Nonlinear interfacial waves in a constant-vorticity planar flow over variable depth 
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Exact Lagrangian in compact form is derived for planar internal waves in a two-fluid system with 
a relatively small density jump (the Boussinesq limit taking place in real oceanic conditions), in 
the presence of a background shear current of constant vorticity, and over arbitrary bottom profile. 
Long-wave asymptotic approximations of higher orders are derived from the exact Hamiltonian 
functional in a remarkably simple way, for two different parametrizations of the interface shape. 
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Large-amplitude internal waves are important phe- 
nomena in near-coastal ocean dynamics (see, e.g., 
Refs.flM^I, and references therein). A nmnber of theo- 
retical models is used to study internal waves analyti- 
cally, including the most popular two- fluid model |9l-[l2|. 
where the system is placed between fixed upper and bot- 
tom boundaries and consists of two fluid layers having 
different constant densities (pi = 1 in the upper layer, 
and p2 = 1 + e in the lower layer). In the simplest vari- 
ant, the rigid walls are both horizontal, and the flow is 
potential within each layer, but additional ingredients 
can be incorporated into the model, such as a back- 
ground shear current and large bathymetric variations 
(see Refs. fTsl - flTj . and references therein). In the general 
case, an analytical study of this system requires quite 
lengthy calculations, especially when fully nonlinear dis- 
persive approximations of higher orders are considered 
[1614241 . In the present work, it will be shown that in 
the case e ^ 1 (the so called Boussinesq limit), there 
exists an elegant and remarkably short way how to de- 
rive fully nonlinear dispersive models of high orders. It 
should be noted that in real oceanic circumstances e is 
indeed small, typically e ^ 0.01, and therefore the theory 
developed below is of practical importance. Moreover, 
our approach easily takes into account the presence of a 
shear background current of constant vorticity, as well 
as arbitrary variations of the bed profile, and that makes 
the model more adequate for studying the real- world pro- 
cesses of interaction of internal waves with tidal currents 
over non-uniform sea bed. We also concern in this work 
the problem of proper parametrization of steep wave pro- 
files and, besides the usual single-valued y = r]{x,t) rep- 
resentation, we consider an alternative parametrization 
which is more suitable to describe arbitrary steep and 
even multi- valued wave profiles. 

Variational formalism. — The method employed here 
is based on a variational formulation for the interface 
dynamics in a two-fluid planar system having constant 
vorticities fli^2 within each layer. Let the (unknown) 
interface proflle be y = rj(x,t), where y is the vertical 
coordinate, x is the horizontal coordinate, and t is the 



time variable. A two-dimensional (2D) velocity field can 
be represented as follows (subscripts 1,2 numbering the 
layers are omitted). 



v(a;, y, t) = nS(x, y) + Vlp{x, y, t), 



(1) 



where the stationary divergence-free field S(a;, y) satis- 
fies equation curl S = 1 and has zero normal compo- 
nent at the rigid boundaries. The unknown potentials 
'Pi,2{x,y,t) satisfy the Laplace equation {dl + dy)(p 1^2 — 
0, with the boundary conditions {'^(pi,2 • n) = at 
a rigid wall (where n the normal unit vector), and 
'Pi,2[x,r]{x,t),t\ = 'ipi,2ix,t) at the interface. Since there 
is a kinematic boundary condition at the interface. 



(vi • n) = (v2 • n) = Vn = m/ V 1 + Vx, 



(2) 



functions ipi {x, t) and ip2{x, i) are related to each other by 
a linear integral operator depending on 77(0;, i). Therefore 
there are only two basic unknown functions, 77(2;, t) and, 
for instance. 



V'(x,i) = ipi{x,t) - (1 + e)i!2{x,t). 



(3) 



It can be proved that the Lagrangian functional for the 
system under consideration has the following structure. 



•^general 



ipVtdx +^ ^-ntdx - H{t], ij}, (4) 



where 7 = {l + e)il2 — ^i, and the Hamiltonian functional 
is the total energy of the system — kinetic plus poten- 
tial. The general idea of the proof is similar to Refs.[25l- 
[27I , and it is based on a generalization of the Bernoulli 
equation for constant-vorticity 2D incompressible fiows 
(subscripts 1,2 are omitted). 



dtf + ne 



vV2- 



gy + P/p = 0, 



(5) 



where Q{x, y, t) is the total stream function in the layer, g 
is the gravity acceleration, and P{x, y, t) is the pressure. 
At the interface, Q and P should be continuous, and this 
requirement gives us the dynamic boundary condition 
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(l + e)9((^2- 



-9t(^i+7e + (l + e)^ 



-egri 



= 0. 

(6) 
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Since, by definition, tlie relation 

dtlllia = [dtVl,2 + ritC)y^l,2\\y=r, (7) 

takes place, and rjt — dx{Q\y=ri)^ Eq.® is equivalent to 

i't+ld-^Vt^SH/dT^, (8) 

while the kinematic boundary condition can be repre- 
sented as 

- r/t = (9) 

Equations ([8]) and (jO)) are equivalent to Eq.Q. Let us 
for example prove that 14, \/ 1 + — —SH/Si/j (consider- 
ation of the dynamic boundary condition is slightly more 
involved but analogous). Since is present in the kinetic 
energy only, the corresponding variation of the Hamilto- 
nian is 

Sn ^ {vi ■\76ipi)dxdy + {1 + e) {'V2 ■V5ip2)dxdy 

JdI Jd2 

VnVlT^[il + f)5^2-6i;i]dx, (10) 



where we have integrated by parts and used the boundary 
conditions and incompressibility. Taking into account the 
definition of '0, we arrive at the required result. 

Exact Hamiltonian theory in Boussinesq limit. — The 
main technical difficulty for application of the Hamilto- 
nian formalism at finite e is the absence of compact ex- 
pression for the kinetic energy of the flow in terms of 
functions rj and ijj. Therefore, below we consider the case 
e ^ 1 and, for simplicity, Vli = Q2 — resulting in 
7 « 0. What is very essential, in this limit ip ~ ipi — 'ip2, 
which means that —tpxix, t)dx is the strength of singular 
vorticity concentrated at the interface. Since the kinetic 
energy of the system in the Boussinesq limit is given by 
integral ICtotai ~ (1/2) / Q^totaidxdy, and 

(9^ + 9^)9 = -ntotai = -n + ^Jxix,t)S[y~r]ix,t)], (11) 

the above observation allows us to express the kinetic 
energy in a closed form through the Green's function of 
the 2D Laplace operator in the domain between the fixed 
boundaries, with zero boundary conditions. An explicit 
expression for the Green's function is known in terms of 
curvilinear conformal coordinates u and v, with w = at 
the upper boundary, w = — 1 at the bottom, {d'^ + dy)v = 
inside the domain, and u{x,y) being a harmonically 
conjugate for v{x, y) [H, [2^. In other words, x + iy = 
z[u + iv), where z{w) is an analytic function of complex 
variable w = u + iv. We can write the Lagrangian in 
terms of the interface profile v — ~q{u,t) and ip{u,t), 



■0 J (u, q)qtdx - H{%Ij, q], 



(12) 



where J{u,q) — \z'{u — iq)\'^ is the Jacobian. The corre- 
sponding equations of motion are 



The Hamiltonian consists of three parts, H = V + JC + S. 
The potential energy is 



r=^ J[lm z{u - iq)]^Re[dui 



iq)]di 



(14) 



/C is the kinetic energy in the absence of external current. 
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In 



sinh[(7r/2)(u2 - Ui - i{q2 + qi))] 



sinh[(7r/2)(M2 - - i(g2 - qi))] 
X'ip[ip2duidu2, (15) 



where qi = q{ui,t), 'i/'i — [dui'{u,t)]\u=ui, and so on. 
The functional S takes into account the presence of a 
background shear current. 



5 = y [Cq — ri6s(u, —q)]^j'du + const. 



(16) 



Here C is a total fiux of the current, and 8s(u,ti) is 
the function which satisfies zero boundary conditions and 
equation (9^ -|- dy)Qs = —1- Usually it is assumed that 
the upper rigid boundary is horizontal at y = 0, and that 
results in ?/(u, —v) = —y{u, v). Then the stream function 
Qs{u,v) can be represented as follows (the idea how to 
calculate Qg is similar to Ref.[30,]), 



Qs{u,v) ^ ^^\^y\u,v) 



sinh kv 



sinh k 



y\u,-i)}, (17) 



where k = —idu is the differential operator, which is 
diagonal in Fourier representation. 

Let us say here that the case 7 7^ 0, though more 
cumbersome, can be considered in a similar manner. 

Long-wave approximations. — Thus, all the terms in 
exact Hamiltonian functional have been explicitly speci- 
fied. However, this exact description cannot be applied 
easily for analytical and numerical studies in view of the 
strongly non-local character of the Hamiltonian. This 
situation with internal waves is in contrast with the ex- 
act description of 2D surface waves of constant vortic- 
ity (soj . where at least numerical implementation is sim- 
ple and efficient with a fast Fourier transform. That is 
why the problem of simplified description for moderately 
steep interfacial waves has attracted much attention in 
last years fl6l - [23 |. Here we suggest a simple procedure 
for the Boussinesq case. We see that the only non-local 
part in the total Hamiltonian is K,. To derive simpli- 
fied quasi-local approximations, we will essentially use a 
Fourier transform of the Green's function. 
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2^ 



In 



sinh[(7r/2)(w - is)] 



sinh[(7r/2)(u — ia)] 



F{k,a,s)e^^^^, 



(18) 



where u — U2 ~ Ui, a = \q2 — qi\, s = 92 + 91, and 



J(w, q)qt = STi/SiJj, 



^J{u,q)il;t^Sn/Sq. (13) 



F(fc,«,.) = ^"^^^(^'°^:;;°f (19) 

smn f£ 

This expression can be easily obtained by the known 
methods of contour integration of analytic functions. It 
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is the simplicity of Eq.([T9]) that allows us to derive the 
mentioned dispersive models of high orders for nonlin- 
ear internal waves. The approximations correspond to 
expansion of F in powers of small quantities ka and ks. 
Therefore, the long-wave assumption actually includes 
three different asymptotic regimes: (i) deep-water the- 
ory, when fc ^ 1 but kq <^ 1, (ii) finite-depth the- 
ory, when q <^ \ and fc ^ 1, (iii) shallow- water the- 
ory, when q 1 and fc ^ 1 [for self-consistency, in this 
regime the bottom variations should be long-scaled as 
well, \x"{u)\/x'[u) <^ 1, where x{u) — z{u + 0«) is a 
purely real function]. Thus, we write 



F = 



2 



-k coth k 



-,2n-l 



2{2n- 1)! 

(a2n _ g2«) 



/c2"-icothfc- 



It is not difficult to understand that each term propor- 
tional to a^^+^fc^™ gives zero contribution to the Hamil- 
tonian. The other terms are easily transformed to finite 
sums of products as t/iJ.?™^ ; thus resulting in the 

quasi-local expansion /C = /Ci + /C2 + /C3 + • • • , where 



/Ci = i / qip'^du, 



(21) 



^2 = / ?MgiV'292[fccothyfc]e''=("^-"i)— duidu2 
2 J 27r 

qi}'[k coth k]{qip')du, (22) 



JC3 



1 



4-3! 
1 

12 



dk 

{q,+q2f^[^P',e''^"'-^'^—duidu2 
Zn 



(23) 



1 



4-4! 



[(gi - 92)" - (-71 + q2fW2 
dk 

x[fc3 cothfc]e*'=("^-"i)— dwidwa 
27r 

{qip')'[kcoihk]{q'^il)'ydu. 



(24) 



Analogously, we derive 
1 



ICa 



IC7 = 



240 



(g^^')'V"' + 5(gV)"(#')' 



40(g^^')"('zV'^" 



du. 



(25) 



1 
360 



2 ■ 7! 



3{q^ijy)"[kcothk]{qxP')" 
+10(g V)" [fc coth k] ((?^ (26) 
(gV')"V^'^+ 7(9'^')"' (#')"' 



+ 21(g^^') (9> ) +35(gV)'"(9''A') du(27) 



It is a simple exercise to calculate the variational deriva- 
tives of the obtained approximate Hamiltonian and sub- 
stitute them into Eqs. ([T3|) . Hopefully, with appropriate 
regularization at short scales, the equations will appear 
to be convenient for future numerical implementation, 
since all the linear operators (the differentiations and 
[fccothfc]) are diagonal in Fourier representation. 

The functionals V and S take a simplified form in the 
case of a slowly varying depth, when the series 

z{u-iq) = x{u)-ix' {u)q-x" {u)'^+ix"' {u)'^ + - ■ ■ (28) 



rapidly converges. Then substitution of Eq. 
Eqs.((ni) and (d?]) gives us 



into 



2 



{x'^x")' 



-x'^x'" 
6 



jdu, (29) 



S « J^'{Cq+^[q\ 



12 

X — 



-q[x 



(x'Y 



q^-q] 



3 J ' 6 
Note that the Jacobian is this case is 

Jiu,q)^x'^ + [ix"f ~x'x"']q^ 



) jdw. (30) 



(31) 



In the shallow-water regime fccothfc « 1 -|- fc'^/3, and 
K. takes a purely local form. 



IC. 



local 



i y <Z(1 - q)ij''du 

{{q^j'y[{l~q)i^']'-[q{l-q)i/W}du 



1 

12 



- -J[q'{l-q)i/n{l-q)i^']'du 

- ^ /[(I - q^qy^riqi^ydu + 0{dt}, (32) 

which is apparently symmetric with respect to change 
q — S> (1 — g). This symmetry is present also in the exact 
expression ([T5)) . since if a ^ a and s (2 — s), then 
F ^ F. 

Alternative parametrization. — Strongly nonlinear in- 
ternal waves are known to have the tendency toward over- 
turning their profiles. The above theory can be general- 
ized to admit more steep, and even multi- valued depen- 
dences y = r]{x,t). To explain the basic idea how to 
manage in that case, we consider here the deep-water 
regime in the absence of current. An arbitrary interface 
shape can be represented in a parametric form, 

x = X{a,t)=a + X{(j,t), y = Y{<T,t)<0, (33) 

with a parameter —00 < a < +00 along the curve. What 
is important, smooth functions X{a,t) and Y{a,t) are 
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able to represent quite steep wave profiles. The corre- 
sponding Lagrangian is 

/Iparam = J [X'Yt ~ XtY']H<T " f ^ Y^X'da 

'{X2-X,f + {Y2 + Y,f 



{X2- Xi)2+(y2-Fi)2 



To simplify the non-local term /C, we again use a Fourier 
transform of the Green's function, but in a slightly dif- 
ferent manner, 



1 

47r 



■In 



{a2-<J^+X2-X^f + {Y2 + Y^f 



(35) 



_{a2-<Ji+X2-XiY + {Y2-Y^f 
= j FiK,X2,Y2,X,,Y^)e^-^--^--^^^, 

with [note the appearance of factor exp{iK(X2 ~ -^1)}] 

-^-\k\\Y2-Yi\ _ p-\K\\Y2+Yl\ 



F 



^ ^iK{X2-Xi) 



2M 



(36) 



Now we expand the exponents in Ea. (l36p in powers of 
the arguments, and take into account that terms propor- 
tional to {X2 - ^i)"|i^2 - yi|2»+ifc"+2" give zero con- 
tribution to the Hamiltonian. As the result, we obtain 



/ I {Y^/y)\k\{Y^')da 



[XYij'ij" + Xi;'{Y'4>')']da 



/[V'"(rV)' + 3(r'^')'(i"V'')V^ 



+- 



{Yx^^y^," - 2{Yx^p'y{xiJ'y 



^{Yi/nx^^'y 



da 



+ / Y4;'\k\{XY4;'yda + 0{d^}, 



(37) 



where \k\ — \da-\ is a pseudo-differential operator. 

It should be noted that if the parametrization of in- 
terface is arbitrary, without any relation between X{a, t) 
and Y{a,t), then the tangential component of the inter- 
face motion (the combination X'Xt +Y'Yt) is not deter- 
mined by variational equations of motion, and therefore it 
remains arbitrary. One can somehow fix the parametriza- 
tion, for example by relation X{a) = —HY{a), where 
H = i sign k is the Hilbert operator. Such parametriza- 
tion is used in the theory of surface waves at the deep 



water (see Refs.[3lH33[, and references therein). For the 
present problem this choice has no special meaning, it 
is only important that arbitrary 2D curves can be rep- 
resented in this way. Then equations of motion for the 
two basic functions Y{a,t) and 'ip{a,t) can be written in 
a non-canonical form (for technical details, see Ref.[33]), 



Yt 



Im 



{^Z'{ 



l + iH) 



{SJC/Si;)/\Z'\']}, 
i/jt Ini{(l-iif) 2{6IC/6Z)Z' + {6IC/6'tlj)i'' 

-eglmZ- ^j'H \{SIC/6^j)/\Z'\^] , (38) 

where = do-ijj{(J,t), Z' = drrZ{<T,t), and 
Z{a, t) = X{(j, t) + iY{(j, t)=(j + {i- H)Y{a, t), (39) 



2{5iC/5Z) ^ (SiC/SX) - i{SiC/SY). 



(40) 



With appropriate regularization at very short scales, the 
above equations can be efficiently simulated on computer 
using fast Fourier transform routines. It should be noted 
that for correct treatment of quite steep waves, say when 
kY ^ A with A ^ 1, one has to expand F in Eg . ((36)) 
up to a sufficiently high order N satisfying the condition 
(2A)^+V(iV + 1)! < 1- For example, the third-order 
approximation Eg. ([57)1 can be good only up to kY « 1. 

Summary and discussion. — To summarize, in this 
work a simple method has been suggested for derivation 
of higher-order dispersive approximations in the theory of 
fully nonlinear planar interfacial waves with a small den- 
sity jump, in a constant-vorticity flow over non-uniform 
bed. Explicit quasi-local expressions have been presented 
up to the 7th order. The method is based on a variational 
formulation of the interface dynamics, and it uses an ex- 
pansion of a Fourier transform of the Green's function 
entering a non-local part of the exact Hamiltonian. 

It should be noted that analogous expansion of a 
Fourier-transformed Green's function can be used also in 
three-dimensional case for the deep-water purely poten- 
tial regime. Another possible application concerns the 
recently suggested two-layer compressible atmospheric 
model [131. Details of the corresponding studies will be 
published elsewhere. 
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